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Correlated imaging, quantum and classical 
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We analytically show that it is possible to perform coherent imaging by using the classical cor- 
relation of two beams obtained by splitting incoherent thermal radiation. The case of such two 
classically correlated beams is treated in parallel with the configuration based on two entangled 
beams produced by parametric down-conversion, and a basic analogy is pointed out. The results 
are compared in a specific numerical example. 



The topic of entangled imaging has attracted noteworthy attention in recent yearsfl H H H H 0,0. This 
tecnique exploits the quantum entaglement of the state generated by parametric down-conversion (PDC), in order 
to retrieve information about an unknown object. In the regime of single photon pair production of PDC, the 
photons of a pair are spatially separated and propagate through two distinct imaging systems. In the path of one 
of the photons an object is located. Information about the spatial distribution of the object is not obtained by 
detection of this photon, but rather by registering the coincidence counts as a function of the other photon position 
0- 0- 0- Gl 01 ■ I n the regime of a large number of photon pairs, this procedure is generalized to the measurement of 
the signal- idler spatial correlation function of intensity fluctuations [6j . Such a two-arm configuration provides more 
flexibility in comparison with standard imaging procedures, as e.g. the possibility of illuminating the object with one 
light frequency and performing a spatially resolved detection in the other arm with a different light frequency, or of 
processing the information from the object by only operating on the imaging system of arm 2 0,0]. In addition, it 
opens the possibility of performing coherent imaging by using, in a sense, spatially incoherent light, since each of the 
two down-converted beams taken separately is described by a thermal-like mixture and only the two-beam state is 
pure(see e.g. [j| and 0). 

In this paper we show that it is possible to implement such a scheme using a truly incoherent light, as the radiation 
produced by a thermal (or quasi-thermal) source. A comparison between thermal and biphoton emission is performed 
in 0, where an underlying duality accompanies the mathematical similarity between the two cases. Here, we consider 




FIG. 1: Correlated imaging with incoherent thermal light. The thermal beam a is splitted into two beams which travel through 
two distinct imaging systems, described by their impulse response functions h\ and /12. Arm 1 includes an object. Detector 
Di is either a point-like detector or a bucket detector. Beam 2 is detected by an array of pixel detectors, v is a vacuum field. 



a different scheme (FigQ, appropriate for correlated imaging, in which a thermal beam is divided by a beam-splitter 
(BS) and the two outgoing beams are handled in the same way as the PDC beams in entangled imaging. A basic 
analogy between the PDC and the thermal case emerges from our analysis. 

Currently there is a debate whether quantum entanglement is necessary to perform correlated imaging [1101. 
The discussion became very lively after the experiment of Q reproduced the results of a ghost image experiment Q 
by using classically correlated beams. We will show here that the spatial correlation of the two beams produced by 
splitting thermal light, although being completely classical, is enough to qualitatively reproduce all the features of 
the entangled imaging. 
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For the sake of comparison we will treat in parallel the cases of entangled beams and of thermal light. For simplicity, 
we consider only spatial variables and ignore the time argument, which corresponds to using a narrow frequency filter. 
We will come back to this point in the final part of the paper. In addition, we assume translational invariance in the 
transverse plane, which amounts to requiring that the cross-section of the source is much larger than the object and 
all the optical elements. In a future publication we will release these assumptions. 

In the entangled case, the signal and idler fields are generated in a type II crystal by a PDC process. Our starting 
point are the input-output relations of the crystal, which in the plane- wave pump approximation read 0, llOl 

h{q) = U i {q)a i {q) + V i {q)a]{-q) i^j = 1,2. (1) 

Here, bi(q) = J ^e~ 1 ^' x bi(x), where bi(x), are the signal (i = 1) and idler (i = 2) field envelope operators at the 
output face of the crystal (distinguished by their orthogonal polarizations), x being the position in the transverse 
plane, , i — 1,2 are the corresponding fields at the input face of the crystal, and are in the vacuum state. The gain 
functions Ui , V% are for example given in [lfj . 

In the thermal case, we start from the input /output relations of a beam splitter 

b\(x) = ra(x) + tv(x) , 62(2?) = ta(x) + rv{x) , (2) 

where t and r are the transmission and reflection coefficients of the mirror, a is a thermal field and v is a vacuum field 
uncorrelated from a. We assume that the thermal state of a is characterized by a Gaussian field statistics, in which 
any correlation function of arbitrary order is expressed via the second order correlation function |l2l | : 

(at(f)o(f')> = J (^ e ~ iq_ ' (2 ~ 2 Hn{q))th (3) 

where (n(q)) t h denotes the expectation value of the photon number in mode q in the thermal state, and we implicitly 
used the hypothesis of translational invariance of the source. In particular, the following factorisation property holds 

(: a\x)a{x')a\x")a{x"') :) = {a\x)a{x')){a\x")a{x"')) + { a\x)a(x' ")){a\x ")a{x') ) , (4) 

where : : indicates normal ordering. A field with these properties is described by a thermal-like density matrix of the 
form 

TT f V- [«q))t h ] m , J 
^ = II 1 L [i + { n(q)) th }^ l TO 'fl< TO 'g1f < (5) 



K m— 



where |m, q) denotes the Fock state with m photons in mode q. 

In both the PDC and thermal case, the two outgoing beams travel through two distinct imaging systems, described 
by their impulse response functions h\{x,x'), h2(x,x') (see Fig. |y. Arm 1 includes an object. Beam 1 is detected 
either by a point-like detector D\, or by a "bucket" detector, which collects all the light in the detection plane, in 
any case giving no information on the object spatial distribution. On the other side, detector D2 spatially resolves 
the light fluctuations, as for example an array of pixel detectors. The fields at the detection planes are given by 

Ci(xi) = dxlhi(xi,Xi')bi(xi') + Li(xi) i = 1,2, (6) 



where L\, Li account for possible losses in the imaging systems, and depend on vacuum field operators uncorrelated 
from 61,62. Since they do not contribute to the normally ordered expectation values that we will calculate in the 
following, their explicit expression is irrelevant. Information about the object is extracted by measuring the spatial 
correlation function of the intensities detected by D\ and D2, as a function of the position a?2 of the pixel of D^,: 

{h{xi)I 2 {x 2 )) = {c\{xi)ci{xi)c\{x2)c 2 {x2)) ■ (7) 

All the object information is concentrated in the correlation function of intensity fluctuations: 

G(x u x 2 ) = {h{x X )h{X2)) - {h{x X )){h{X2)) , (8) 
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where (J^(afj)) = {c\{xi)ci(xi)) is the mean intensity of the i-th beam. Since c\ and cj commute, all the terms in Eqs. 
0,© are normally ordered and Li,L2 can be neglected, thus obtaining 



G(xi,X2) = / dx{ j Ax" J dx' 2 J dxghKxi, x")hi(xi, x^h^i^i, x%)h2{x2, x'2) 
\bl(x[')b 1 (x{)bl(x^b 2 (x^)) - (blixl'Mxl)) (bUxZMx^)) . 



(9) 



In the thermal case, by taking into account the transformation J5J and that v is in the vacuum state, 61 and &2 in 
Eq.@ can be simply replaced by ra and ta. Next, by using Eq.Q, we arrive at the final result 



G(xi,X2) = \tr\ 



Ax[ I Ax2h\{xi,x[)h 2 {x2,X2)(a\x[)a(x2)) 



(10) 



On the other hand, also in the PDC case the four-point correlation function in Eq.@ has special factorization 
properties. As it can be obtained from Eq.Q [1(1 ]. 



(bUx^hixDbUx^'M^)) = (b\&'M2l)){bU22)h&)) 

+ (b\W)b\{^))(b x {x*[)b2&)) 

By inserting this result in Eq. @, one obtains 

dx[ I Ax2h\{xi,S[)h2{x2,x^{bx{x[)b2{x2)) 



G{xi,x 2 ) - 
where, by using relations l|T]l. 



dq 

(2nf 



(11) 



(12) 



(13) 



There is a clear analogy between the results in the two cases. Apart from the numerical factor \tr\ 2 and the presence 
of h\ instead of hi, the second order correlation (a^ (x)a(x')) and the function {n(q)) t h play in Ea. l|10[l the same role 
as the correlation function (b\(x)b2{x')) and U\{q) Vi(—Q) in Ea. (|12f> . Most importantly, in both Eas-ljlOfland l|12() the 
modulus is outside the integral, a feature that ensures the possibility of coherent imaging via correlation function (se e 
e.g. The correlation function (a) (x)a(x')) governs the properties of spatial coherence of the thermal source [l^|. 
The correlation length, or transverse coherence length l co h, is determined by the inverse of the bandwidth Aq of the 
function (n(q)) t h- The same holds for the correlation (bi(x)b2(x')), and the function Ui(q)V 2 (—q) in the entangled 
case. 

Let us now analyse two paradigmatic examples of imaging systems, mutated from the discussion of 6], and de- 
scribed in FigEl In both examples the set-up of arm 1 is fixed, and consists of an object, described by a complex 




FIG. 2: Imaging scheme. L denotes two identical lenses of focal length /. The distance z is either z — f or z = 2f. 

transmission function T(x), and a lens located at a focal distance / from the object and from the detection plane. 
Hence, hi(xi,x{) — — jj exp ^— Tf^i ' &i)T(x[), with A being the wavelength. In arm 2 there is a single lens placed 



4 



at a distance z both from the source and from the detection plane 2. 

In the first example we assume z=f (we take the two lenses identical for simplicity), so that h 2 (x 2l X2) = 

— jj exp ^— ^fx 2 • x£\ ■ By inserting these propagators into Ea. (jTU|) . and taking into account the expression of 
(a* (x)a(x')) given by Eq.JSJ, we obtain 



G(x u x 2 ) = 



(A/) 4 



(14) 



where T(g) = / ^e~ 1 ^' yi T{x) is the amplitude of the diffraction pattern from the object. This has to be compared with 
the result of the entangled case (see Eq.7 of |(J), where the combination x 2 + x\ appears instead of x 2 — x% and XJ\V% 
instead of (n) t h- In both cases, the diffraction pattern from the object can be reconstructed, provided that the spatial 
bandwidth Aq is larger than the maximal q vector appearing in the diffraction pattern, or, equivalently, provided 
that l co h < l , where l is the smallest scale of variation of the object spatial distribution. Best performances of the 
scheme are achieved for spatially incoherent light, l co h — + 0. We remind the well-known result that, when l co h < l , 
no information about the diffraction pattern of the object is obtained by detection of light intensity distribution in 
arm 1 by an array of pixels. 

In the second example, we set z=2f, so that h 2 (x 2 , x 2 ) = 8{x 2 + S 2 ) exp ( — i|x2 1 2 XT ) - Inserting this in Eq. I|10|) . we 
get: 



G(xi,x 2 ) 



\rt\ 



(A/) 2 
\rt\ 2 
(A/) 2 



[n(xi^)) th 



\T(-x 2 ) 



(15) 
(16) 



where in the second line l co h < l was assumed, so that (n(q))th is roughly constant in the region of q plane where the 
diffraction pattern does not vanish, and it can be taken out from the integral in i|15|) . In this example the correlation 
function provides information about the image of the object. A similar result holds for the case of entangled beams 
(see Eq.(8) of @). 

Our results appear suprising, if one has in mind the case of a coherent beam impinging on a beam splitter, where 
the two outgoing fields are uncorrelated, i.e. G{x\,x 2 ) = 0. However, when the input field is an intense thermal 
beam, i.e. the photon number per mode is not too small, the two outgoing field are well correlated in space. To prove 
this point, let us consider the number of photons detected in two small identical portions R ("pixels" ) of the beams in 
the near field immediately after the beam splitter, TV; = J R dxbl(x)bi(x) i = 1,2, and the difference 7V_ = Ni — N 2 . 
Making use of the transformation (J2J, it can be proven that, for \r\ 2 = \t\ 2 = 1/2, the variance (5N^_) = (iV 2 ) — (TV_) 2 
is given by 



(SN 2 _) = (N,) + (N 2 ) , 



(17) 



which corresponds exactly to the shot noise level. On the other side, by using the identity (<5iV 2 ) = (5N 2 ) + (SN 2 ) 
2(5Ni8N 2 ) and taking into account that (5N 2 ) = (5N 2 ) for |r| 2 = \t\ 2 , the normalized correlation is given by 



C 



dcf (SNJNz) = (/Vi) 



(18) 



For any state < \C\ < 1, where the lower bound corresponds to the coherent state level, and the upper bound is 
imposed by Cauchy-Schwarz inequality. For the thermal state, provided that the pixel size is not too large with respect 
to l coh , (8Nf) « (Ni) + (iVi) 2 , so that the correlation (TSJ} never vanishes. For thermal systems with a large number 
of photons (N\ ) I {5N 2 ) w l/(Ni) << 1, and C can be made close its maximum value. Even more important, it is not 
difficult to show that Eas. l|17fl and (|18fl hold in any plane linked to the near field plane by a Fresnel transformation, 
in the absence of losses. In particular, a high level of pixcl-by-pixcl correlation can be observed in the far-field plane. 
We remark that, despite C can be made close to 1 by increasing the mean number of photons, the correlation never 
reaches the quantum level, as shown by Ea. (|17|) . 

For the entangled beams produced by PDC, spatial correlation is present both in the near and in the far field, with 
the ideal result (<5iV 2 ) = 0, C = 1 in both planes In this case, the far field correlation is between symmetric 

pixels, due to the different way in which the four-point correlation function factorises. In jfj we analysed the effect 
of replacing the pure PDC entangled state with two mixtures that exactly preserve the spatial signal-idler quantum 
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correlations, either in the far or in the near field. It turned out that, when considering the "far-field mixture", the 
pure state results in the z=f configuration of Fig[2]can be reproduced, but there is no information about the image 
in the z=2f configuration. The converse is true considering the "near-field mixture" . The two beams generated by 
splitting thermal light are instead imperfectly correlated both in the near and in the far-field. However, by using 
intense thermal light, classical intensity correlation is strong enough to reproduce qualitatively the results of both the 
z=f and the z=2f configuration. 

A complete comparison of the performances in the classical and quantum regimes requires extended numerical inves- 
tigations, describing realistic thermal sources, which is outside the scope of this paper. A key role in this comparison 
will be played by the issue of the visibility of the information, which is retrieved by subtracting from the measured 
correlation function {JJ, the background term (Ii(xi))(l2(x2)} (see Eq.J^J). A measure of the visibility is given by 

1 ; _ G(S 1 ,x 2 ) 

V - </l(?l)/2(? 2 ))- 

A first remark concerns the presence of (n(q))th in Eq. 1)10(1 in place of U\(q)V2{—q) in Ea. ((12|l . As a consequence, in 
the thermal case G(x\, X2) scales as (n(q)) 2 h , while in the entangled case, it scales as \Ui(q)V2(— q)\ = (n(q)} + (n(q)) 2 , 
where (n(q)) = \V 2 (-q)\ 2 = | Vi (g)| 2 is the mean number of photons per mode in the PDC beams, and |[/i(<f)| — 
l + \Vi(q)\ 2 (see e.g 10]). This difference is immaterial when (n(qj) >> 1, while it becomes relevant for a small photon 
number, because the background (Ii(xi)} {hi^)} on (n(q)) 2 is negligible with respect to G{xi,x~2) oc (n(q)). Hence, 
in the regime of single photon pair detection, the entangled case presents a much better visibility of the information 
with respect to classically correlated thermal beams (see also (lif). 

A second remark concerns the role of the temporal argument. Calculations that will be reported elsewhere show 
that the visibility scales as the ratio between the coherence time of the source T co h and the detection time (see also 
Ilflj). In practical applications, this implies that conventional thermal sources, with very small coherence times, 
are not suitable for the schemes studied here. A suitable source should present a relatively long coherence time, as for 
example a sodium lamp, or the chaotic light produced by scattering a laser beam through a random medium UM. As 
a special example of a thermal source, one can consider the signal field or the idler field generated by PDC. Fig. 



(a) Entangled 

(b) Thermal 




x 2 /x Q 

FIG. 3: Numerical simulation of the reconstruction of the diffraction pattern of a double slit in the scheme z=f of Figl2l 
G(xi,X2) versus £2, after 10000 shots for (a) entangled signal/idler beams from PDC, (b) classically correlated beams by 
splitting the idler beam, (c) is the analytical result of Ea, 1)14(1 . Parameters are those of a 4 mm /3-barium-borate crystal 
(Icoh = 16.6/im, T co h = 0.97ps). The pump waist is 664/xm, and the pulse duration is 1.5ps. xq = AqXf /(2tt) 

13 shows the results of a numerical simulation for the reconstruction of the diffraction pattern of a double slit, in the 
scheme z=f of Fig|3] It compares the use of signal/idler entangled beams (curve (a)), and two classically correlated 
beams obtained by symmetrically splitting the idler beam (curve (b)). Parametric gain is such that (n(q)) « 750 at 
its maximum for (a), and (n(q)) ~ 1500 for (b), so that the the mean photon number of beams 61 62 is the same in the 
two simulations. This example not only shows that in the regime of high photon number the quantum and classical 
correlation offer similar performances, but also that although the visibility is rather poor, V « 0.05 in both cases, this 
is enough to correctly retrieve the desired information after a reasonable number of pump shots. The numerical data 
for the image reconstruction (scheme z=2f of Fig|5J confirm these results. 

In conclusion, we have shown a deep analogy between the use of entangled beams and classically correlated beams 
from a thermal source in imaging schemes based on correlation measurements. As it was already recognized in other 
contexts (see e.g. 16] ), in the small photon number regime, a definite advantage of the quantum configuration is 
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represented by a better visibility. 

Our result, that is possible to perform coherent imaging without spatial coherence by using thermal light in combi- 
nation of a beam splitter is reminiscent of the Hanbury-Brown and Twiss interferometric method for determining the 
stellar diameter |l7|. or of the detection of the fringes arising from interference of two independent thermal sources 
(lif. However, here, we define a technique to achieve a full coherent imaging with a great deal of flexibility. Since the 
required correlation is classical, a high quantum efficiency of detectors is not necessary. 

This work was carried out in the framework of the FET project QUANTIM of the EU. We are grateful to S. 
Sergienko, B. Boyd and E.Lantz for stimulating discussions. 
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